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SUMMARY 

? H — 

The components of the planetary perturbations in rectangular co- 
ordinates are given as series expansions in the well-known Hansen 
components. The results provide formulas for transforming a general 
theory in one set of components into the corresponding general theory 
in the other set of components by simple series multiplications. The 
results suggest a useful method of improving the convergence of suc- 
cessive approximations when computing perturbations in the rectangu- 
lar coordinates. 
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Musen’s method of planetary perturbations (Reference 1) bears the greatest resemblance to 
the method of Hansen (Reference 2). The relations between the perturbations computed from the 
two methods have been studied extensively by the author, since Hansen’s method has been used in 
many cases. A comparison of the integrands shows that the first-order results should be identi- 
cal, term by term. There is no need to give this comparison here, however, since the main 
interest is in the higher order differences. These higher order results reveal the small differ- 
ences in the coefficients and suggest a modification which has been used to improve the conver- 
gence of successive approximations in Musen’s method. 

In Hansen's form, the rectangular coordinates referred to the mean plane of the orbit at the 
fundamental epoch are given by 


X = r COS (f + cr) cos b , 

y = r sin (f + cr) cos b , 

z = r sin b , 

using the notation of Clemence (Reference 3, page 317). Here, cr represents the reduction of the 
orbital longitude from the osculating plane to the mean plane. The approximation 

_ P 1 ds 

- r ■ T s d7 

has been used by Hill (Reference 4) and Clemence (Reference 3). The other quantities are defined 
in terms of the elliptic elements and the perturbation components by the following set of formulas: 

cr — n + ri ^ t — t ^ 

D t>0 l 0 ) * 

*In Part I of this report, a description of a program for computing Hansen’s planetary perturbations was given. Examples were included 
for six minor planets. 
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e - e sin e g + nSz , 

T = a( 1 * e cos <f) , 
r = T (1 + v) , 

T cos f = a (cos T * e) , 

T sin f = a y 1 - e 2 sin T , 

a 

sin b = *=■ u . 


Using these formulas, 

x = a( 1 + v) j\cos e - e) cos a - yi - e 2 sin T sin a 
y = a(l + v) [/i - e 2 sin T cos a + (cos F - e) sin crj 


- z(l + v) u . 


Writing Musen's formula for the perturbations in the form 


1 dr o 


6r - at Q + 0 [— dt I + 7 (aR) 


where 


r Q = aP(cos e - e) + a y 1 - e 2 Q sin e 


the corresponding components are 


o a 

X - a( 1 + a )(cos e - e) - a/3 r sin € , 

o 


= a(l + a) fT- e* sin e + a/3 — - e 2 cos e 


ay 


From the latitude components, 


y - ( 1 + v) u 


without approximation. 


] cos b t 
cos b , 
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Putting 


Se - 6-e and r Q - a(l - e cos e) , 


the remaining components yield 


1+a - (1 + v) cos b <cos a + cos a — (cos 8e - 1) 

' o 


~ sin a j/l - e 2 — j^sin $e + — sin e (1 - e cos e) ^ 


and 


J3 - (1 + v) cos b cos a [(1 - e cos e) sin Se + e sin e (1 - cos Se)] 

e e "|| 

cos 8e t (1 - e cos e) cos T r (cos € - e) > 

1 ~ e* 1 “ e JJ 


+ sin cr y 1 - e 2 


These formulas are exact. The remaining steps are to expand cos b in powers of u, cos cr and 
sin a in powers of a, and 8e in powers of nSz. In performing these expansions, only those terms 
which appear to be of significance will be kept. The results can easily be extended to include 
higher powers of the perturbations if necessary. The largest component is generally nSz, and the 
fifth power will be retained. The next largest component is v , which appears here as a linear 
factor, so that no expansion in powers of v is necessary. The third component, u , is generally 
much smaller. It enters through the factor cos b , so that only the even powers of u occur, u 2 is 
kept here; but u 2 nSz, u 2 v, and u 4 are dropped. Finally, a- is taken to be quite small so that cr 2 , 
a • n8z, and a * v are dropped. The fact that terms of the fifth order in nSz may be important 
is illustrated in the case of the minor planet (1373) Cincinnati (Reference 5). This degree of 
approximation should be adequate for many cases of interest. Now, 


cos b 



cos cr - 1 , 


sin cr - cr 


and 


1+a 


(1 + v) 


1 + — (cos 8e 
o 


~ 1) 



sin 6 
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P = (1 + v) sin 8e + e sin e (1 - cos Se)J + , _ =$ J 


From Kepler's equation 


8c - nSz + e (sin £ - sin £) 


which gives 


— nSz + 
1 o 


- T sin € nSz 2 + |-y sin 


6 n8z 3 


5e 3 / a \ 7 , 5e 2 / a \ 6 . e / a\ 5 . 4 

— * — I — i cin 3 f + - 4 n I— — ! cos e sin e + 777F - ! sin £ nSz 


"T (77J sin e + ~TT cos 6 sin e + I? \ r 0 


7e 4 / a 


/_ sin 4 e -- g - COS 6 sin € + 12 ^ 


4 (4) sin2 £ + T5o (4) 


— ] cos € n8z 5 + 0 (n8z 6 ) 

r o / 


Using the resulting expansions of cos Se and sin Se , 


a - v + (1 + v) J- -J 4) "I* 2 + t(t;) sin e nSz 3 


■^fey - 41 fey * m\ 


+ [- 7e(1 s ' e -- (4) sin e + 4 ( 4 ) sin e (4) sin e 


e . 1 ( a \ 2 2 

W ^ cr sin £ “ -o* T u 









Next nSz, v, and u are given as expansions in a, 0 , and y and the small quantity a by inverting 
these expansions. The results are 


V 



“V'T 77 


(1 


2a + 3a 2 ) fi 3 


+ 



(1 - 3a) /3 4 


+ 


7e( 1 - e 2 ) 
8 



e 



+ 


e ( a \ 7 

~ sin e 
Z \ 0/ 


/3 s 


f' 


- e 2 


sm 6 - 


nSz = (l 


a + a 2 


+ a 4 )/ 3 * T fe-) 


(l - 3a + 6 a 2 ) yS 3 + * 5 " (— J sine(l“4a)y6' 


+ 







u 


yS 2 + a 2 


Terms of the third and higher orders in a and yS factored by y in the last equation have been 
dropped. 

If the perturbations in Hansen f s components are computed from mean elements, then certain 
terms in nSz and u are zero. The corresponding terms in fi and y are of the second order. 

It was mentioned before that the first-order results from the two methods are identical and 
therefore may be taken to represent either a , /3 , y or v, nSz , u. In either case, 


a 


v 



fi - nSz - v * n 8 z + • • • , 


y - u - v * u , 


where the omitted terms are of the third order. When nSz (or /3) contains large terms resulting 
from near resonance, there will be a significant difference between a and v in the constant terms 
and the terms with twice the long-period argument. The importance of this difference comes 
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from the strong cross-action of these terms in the radius vector with the long-period terms in 
the longitude in the higher approximations. The convergence of successive approximations in 
Musen T s method is considerably improved if the appropriate changes are made to these terms in 
a . If these changes are not made, they will appear in the next approximation anyway; but the new 
values of the long-period terms in fi will be quite bad. 

These considerations with respect to the long-period terms also should be taken into account 
in the application of other methods using rectangular coordinates. 

The author expresses his appreciation to Dr. Peter Musen. who suggested the work on general 
perturbations, and to Neal Paris for his assistance in programming. 


(Manuscript received February 18, 1965) 
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